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Abstract 

^  A  variational  approach  is  outlined  for  the  deduction  of  impurity 
confinement  scaling  laws.  Given  the  forms  of  the  diffusive  and  convective 

components  to  the  impurity  particle  flux7"~-we^  present;  a  variational 

'A 

principle  for  the  impurity  confinement  time  in  terms  of  the  diffusion 
time-scale  and  the  convection  parameter,  which  is  a  non-dimensional 
measure  of  the  size  of  the  convective  flux  relative  to  the  diffusive  flux. 
These  results  are  very  general  and  apply  irrespective  of  whether  the 
transport  fluxes  are  of  theoretical  or  empirical  origin.  The  impurity 
confinement  time  scales  exponentially  with  the  convection  parameter  in 
cases  of  practical  interest.  -■  ‘  “ 
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1 .  INTRODUCTION 


Experimental  studies  of  trace  impurity  transport  (HAWKES  et  al. 

(1987),  MARMAR  et  al.  (1982)  represent  tvpical  examples)  provide  a  useful 
way  of  probing  tokamak  plasmas.  The  present  article  is  concerned  with 
certain  general  features  of  impurity  confinement  scaling.  We  consider  a 
cylindrical  model  of  tokamak  magnetic  flux  surfaces  for  simplicity.  The 
trace  impurity  number  density  n^(r,t)  is  assumed  to  satisfy  the  following 
transport  equation. 

3nT  1  A 

— i  -  -  1  _  (r  Tjlr.t))  (1) 

dt  r  dr 

In  practice,  the  impurities  can  exist  in  various  ionization  states  with 
different  values  of  .  In  principle,  one  must  then  solve  a  series  of 
coupled  equations  of  the  type  of  Eq.(1)  with  sources  to  model  the 
transitions.  The  theory  presented  in  this  paper  is  applicable  to  this 
more  general  case  also.  However,  it  is  best  understood  in  terms  of  the 
simple  model  provided  by  Eq.(1). 

The  impurity  radial  flux  T  (r,t)  is  known  (both  empirically  and 
from  neo-classical  theory)  to  take  the  general  form, 

an 

ryr.t)  =  -  Dj.tr) -  -  Ujtrlnj  (2) 

Sr 

The  diffusivity  D^fr)  and  the  (nominally)  inward  convection  velocity 


-  1  - 
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The  fiiSt  two  conditions  specify  the  nature  of  the  transport 


co-efficients.  It  is  usually  the  case  that  the  profile  functions 
f^(r/a),  f^(r/a)  are  monotonic,  increasing.  Conditions  (3)  and  (4)  are 

the  boundary  and  initial  conditions  respectively.  For  experimental 
applications,  it  is  useful  to  regard  Finxtial^r*  as  a  ^unctaon 


concentrated  near  r  =  a . 


Given  the  above  data,  the  initial  value  problem  for  Eq.(1)  can  be 
solved  numerically  in  general.  In  special  cases,  analytic  solutions  may 
also  be  constructed  in  terms  of  known  functions.  To  discuss  the 
properties  of  n^(r,t)  we  introduce  two  parameters. 


Let 


Tdiff  * 


a2 

DI(a) 


(3) 


a  0  (a) 

and  C  5  - i -  (4) 

3  0  (a) 


Clearly  is  a  “easure  of  the  impurity  confinement  time  (to  be 

defined  precisely)  in  the  absence  of  convection,  while  C^  is  a 
dimensionless  measure  of  the  convective  flux  relative  to  the  diffusion 
flux.  C^  is  usually  referred  to  as  the  Peclet  number  in  fluid  mechanics 
(JERRARD  and  McNEILL,  1986)  though  it  is  convenient  for  our  purposes  to 
call  it  the  convection  parameter.  From  n^(r,t)  we  may  form  the  line 
average 


N(t) 


I  n  <r,t)  — 

0  a 


(5) 
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A  plot  of  this  iFig.  1 )  shows  a  vary  sharp  rise  followed  by  a  nearly 
exponential  decay  in  time.  The  experimental  brightness  function  is 
related  to  i3(t),  though  not  simply.  However,  the  decay  time  is  a 

measure  of  the  impurity  confinement  time. 

From  dimensional  analysis,  it  is  obvious  that  the  following  "scaling 
law"  must  hold. 


i. 

imp 


/t 


diff 


F(Ca) 


(6) 


Fig.  2  shows  the  results  of  a  set  of  numerical  solutions  of  Eq.(1)  for 

typical  forms  of  f  ,  f  ,  F  ...  over  a  range  of  C  .  F(C  )  can  be 
DU  initial  a  a 

written  (for  C^  >  0)  FfC^)  -  ACa  exp[BC^}.  The  constants  A,  B,  a,  p 

depend  only  on  fp,  f^  and  not  on  It  *s  noteworthy  that 

F ( C  )  is  exponentially  dependent  on  C  whereas  Eq.(1)  has 

3i  cl 

co-efficients  Oj  which  only  depend  on  linearly.  This  is  a  direct 

result  of  the  cancellation  between  the  diffusive  and  inward  convective 
fluxes. 


3.  ANALYTIC  THEORY  AND  VARIATIONAL  PRINCIPLE  FOR  THE  SCALING  FUNCTION 


The  above  results  can  be  understood  quite  simply  by  the  application 
of  the  following  theorem. 

Theorem: 

Consider  the  eigen-value  problem  defined  by 


4 


(7) 


—  — (r  D  (r)—  +  r  U  (r)$]  *  -H 
r  dr  dr 


together  with  the  boundary  conditions,  —  =  0  at  r  =  0,  $  =  0  at  r 

dr 


a.  The  eigen-values  {x^  form  a  real,  positive,  monotonically 

increasing  unbounded  sequence,  ie,  0  <  X^  <  X2  <  ■ •  •  <  X  .... 


b.  The  corresponding  eigenfunctions  4>^  may  be  written  in  the  form 

$  “  p( r )  V  (r) 

n  n 

r  U  (r) 

where  p(r)  5  exp  {  -  /  -  dr}  and  V  satisfy 

0  Dr(r)  " 


I  — 

r  dr 


(r  D  p 


dV 

n . 

dr 


X  pV 
n  n 


(8) 


and  are  normalised 

a 

Clearly,  6  =  / 

ran  ' 


a 

by  1  =  /  p(r)V  2  rdr. 
0  n 

p(r)V  V  rdr  holds, 
n  m 


c.  The  solution  to  Eq.(1)  is  given  by  the  eigenvalue  expansion 


“  -  X  t 

n  (r,t)  =  I  e  n  A  $  (r)  (9) 

I  n  n 

n=  1 


where  F  (r) 

initial 


l 

n=1 


A  $  (r) 
n  n 


(10) 


and  t  ~  -1_  . 
imp 

M 


e.  The  eigen-value  is  given  by  the  Rayleigh-Ritz  principle 


V(r)  a 


J  p(r)  V2  rdr 
0 


Where  the  minimum  is  taken  over  real,  continuously  differentiable 
dV 

functions  s.t  —  «  0,  r  =  0;  V  =  0,  r  =  a 
dr 


Proof : 


r  Ut 

The  substitution  <S  =  p(r)V  with  p(r)  =  exp  {-  /  —  dr}  converts 

0  D 

1 

(7)  to  the  self-adjoint  eiqen-value  problem  (the  boundary  conditions  are 
unaffected) 


1  — ( r  DJp  fl)  «  -  \  p  V  (13) 

r  dr  dr 


—  “0  at  r  =  0  and  V  *  0  at  r  =  a  .  The  statements  a,  b,  c,  d 
dr 


6 


and  e  follow  from  the  standard  theory  given  in  the  well-known  texts  of 
CODDINGTON  and  LEVINSON  (1955)  and  COURANT  and  HILBERT  (1953). 

We  now  apply  the  variational  principle  to  deduce  FIC^)  for  an 

experimentally  interesting  case.  The  point  of  this  application  is  this: 

if  we  are  only  interested  in  F(C  ),  and  not  in  n  (r,t),  we  need  not 

a  I 

solve  either  Eq.(1)  or  the  eigenvalue  problem.  We  simply  take  a  suitable 
trial  function  V  and  evaluate  X^  approximately  by  calculating  the 
Rayleigh-quotient 


7 


Putting  V  =  (1-x2)  and  evaluating  the  Rayleigh-quotient  we  get, 


F ( 0 )  *  2  *  0.167 
6 


In  the  general  case,  we  consider  the  trial  function  (other  choices  are 


possible),  V  =  (l  -  exp  {_£  (x2-1)j  which  satisfies  the  boundary 

C  2 

a 

conditions  and  reduces  to  the  previous  one  in  the  limit  C  -►  0 .  The 
Rayleigh  quotient  is  readily  evaluated  in  closed  form.  We  then  obtain  the 
result 


x. 

imp 

Xdiff 


where, 


R(x) 


4xe 


-2x 


(U-1)e*  + 
{1  -  2xe'* 


For  large  C , 


(  14) 


In  Fig. 3,  we  have  plotted  the  points  obtained  by  a  fine-grid  numerical 
solution  of  Eq.(1)  and  the  curve  given  by  Eq.{14)  for  logt^^FtC  ). 


3 


4. 


DISCUSSION  AND  CONCLUSIONS: 


The  above  results  demonstrate  that  for  given  profiles  f  and  f  , 

DU 

the  confinement  time  t,  is  determined  by  two  parameters.  These  are 
imp 


the  diffusion  time-scale  t. 


—  and  the  convection  parameter. 


a  U  (a) 

C  =  _  .  In  a  sense,  diffusion  'always  wins'  no  matter  how  large 

3  Djla) 

C  is.  This  is  due  to  the  fact  that  A,  >  0  however  large  C  is. 


Clearly  this  is  an  effect  of  a  finite  domain  and  the  condition 


n^(a,t)  =  0.  However,  if  C  >  5  say,  the  exponential  dependence  on 

C  asserts  itself  and  z.  can  be  very  much  larger  than  The 

a  imp  diff 

calculations  also  show  that  an  increase  or  decrease  of  z.  can  be 

imp 

mediated  solely  by  the  variation  of  C  without  t. ...  tie  D  )  itself 

a  diff  I 

changing. 


Experimental  results  on  trace  impurity  transport  can  be  simply 

interpreted  in  terms  of  the  above  model.  In  Alcator  C  for  example, 

assuming  D  to  be  anomalous  (of  order  y  )  and  U  to  be  the 

Pf irsch-Schluter  convection,  we  (HAAS  and  THYAGARAJA,  1987)  have  been  able 

to  explain  many  features  of  the  results  obtained  by  Marmar  et  al .  In 

contrast,  DITE  (AXON  et  al,  1987;  HAWKES  et  al,  1987)  can  operate  in 

apparently  two  regimes  depending  on  the  plasma  density  rate  of  increase. 

This  could  be  pictured  as  being  due  to  the  variation  of  C  .  Pursuing 

a 

this  idea  further,  our  phenomenological  model  predicted  a  strong 
degradation  of  impurity  confinement  due  to  auxiliary  heating.  This  effect 
is  a  consequence  of  the  general  form  of  the  scaling  function  F  and  the 
decrease  of  suggested  by  the  physics  of  the  model.  The  methods  of 


this  paper  may  be  applicable  more  generally  to  plasma  confinement  and  to 
impurity  confinement  with  sources. 

The  form  of  the  scaling  function  F  cannot  be  deduced  without 

reference  to  the  boundary  value  problem.  In  particular#  no  scaling  or 

group  theoretic  analysis  of  Eq.(1)  can  ever  lead  to  the  actual  form  taken 

by  F  .  This  example  indicates  that  scaling  arguments  based  on  invariance 

principles  alone  should  be  treated  with  caution.  However,  it  is  known 

from  standard  eigenvalue  theory  that  since  the  co-efficients  of  Eq.(3)  are 

analytic  functions  of  C^  (through  p  ),  this  must  also  be  the  case  for 

the  eigen-values .  Thus,  FiC^)  must  be  an  analytic  (in  general 

transcendental)  function  of  C  .  This  is  in  contrast  to  nonlinear 

a 

problems  (eg  plasma  heat  transport)  where  FIC^)  can  depend 

discontinuously  on  C  (i.e.,  exhibit  bifurcations).  Of  course,  x 

a  diff 

and  C  themselves  could  change  discontinuously  in  an  experiment, 
a 

Figure  Captions : 

Fig.  1.  Typical  calculated  variation  of  N(t)  with  time  using  Eq.(1) 
and  Eq. ( 5 ) . 

Fig.  2.  Calculated  variation  of  x.  as  a  function  of  C  for  a 

imp  diff  a 

specific  choice  of  f  and  f  . 

DU 

Fig.  3.  Comparison  between  the  results  of  numerical  solutions  of  Eq.(1) 
(shown  by  open  circles)  and  the  variational  formula  Eq.(14) 

(solid  curve).  Log10  is  plotted  against  C^  • 

xdif  f 
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Fig.  1  Typical  calculated  variation  of  N( t)  with  time  using  Hq.(l) 
and  Eq.(5). 
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Fig.  2  Calculated  variation  of  as  a  function  of  C„  for  a 

specific  choice  of  So  and  A  ■ 
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Fig.  3  Comparison  between  the  results  of  numerical  solutions  of 
Hq.(l)  (shown  by  open  circles)  and  the  variational  formula  Eq.(]4) 

(solid  curve).  Log,0 is  plotted  against  C„. 
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